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We show that ofF-diagonal nearest neighbor disorder in quasi-one-dimensional single particle tight- 
binding coupled chains leads to anomalies in the density of states and in the mean conductance, that 
can be interpreted as due to specific resonances of the band structure of the perfect system underlying 
the disordered one. We demonstrate that this phenomenology may appear not necessarily at the 
energy i5 = as reported so far in the literature and we show that also the even-odd chain number 
effect on the mean conductance is ruled by the same resonances. For different cases we provide a 
series expansion for the resonant contribution to the mean conductance. These expansions allow to 
make predictions well verified numerically. 



The existence of a divergent density of states at the 
center of the spectrum {E = 0) of a one-dimensional 
phonon model with off-diagonal disorder goes back to 
the pioneer paper of Dyson in 1953 Q. The effect of 
this kind of disorder on the divergence of the localization 
length at the energy E = was then evidenced 0, Q , 
in contrast with the general wisdom that in strictly one 
dimensional systems with uncorrelated diagonal or off- 
diagonal disorder all the states are localized 0, Q . 

Recently further interest has grown on off-diagonal dis- 
order in quasi-one-dimensional systems composed of N 
coupled chains. In this case it has been shown that for 
odd N the localization length and the density of states 
diverge at the center of the band, at zero energy, while 
they remain finite for even iv|^. 

It has been claimed that the origin of delocaUzation 
of the states at i? = 0, is the existence of an additional 
(chiral) sublattice symmetry present in system with off- 
diagonal disorder only Q- 

To our knowledge only the point E = has been in- 
vestigated in the literature. In this paper we look for the 
possibility of occurrence of the anomalies described above 
at other energy points. We show that such anomalies can 
arise in every point of the energy spectrum, provided a 
well defined condition on the energy bands of the perfect 
underlying chain system is fulfilled. 

We start considering an infinitely extended periodic 
quantum wire composed by N interacting chains; diag- 
onal and off-diagonal disorder is then introduced in a 
strand of length C of the wire (the lattice constant is 
taken as unit of length) . 

The Schrodingcr equation for the multichain system 
considered is: 



Un'ipn + Tn-li>n-l + T^V'n+l = Elpn, \fn, 



(1) 



where [/„ is the matrix connecting the sites belonging 
to the n-th column and r„ is the matrix specifying the 
hopping amplitudes from the column n to the next one. 



The order of these matrices is equal to the number of 
coupled chains, N. The index n g Z labels the columns 
along the wire and the ipn arc iV-dimensional vectors. 

In order to introduce the disorder we write the matrices 
Un and T„ in the following form: 



Un 



(2) 



where U and T are the matrices, independent of n, defin- 
ing the periodic chain system and and are ran- 
domly distributed matrices. We take ^ only for 
n = and VFj 7^ for n = 1, 1. The dis- 

tinction ^ defines the ordered part and the perturbation 
part (disorder) of the Hamiltonian, namely H = Ho-\-W. 
All the elements of the random matrices have the form 



„i-a, ) 'Xi 



(3) 



where all the x and x are independent and uniformly 
distributed real random variables, with zero mean and 
second moment equal to one. The terms wfj"^^ and cx^j^'^ 
do not depend on the column index n and regulate the 
modulus and the phase of the disorder. 

Let us consider first the periodic chain system de- 
scribed by the Hamiltonian Hq. Let Z{9b) be the matrix 
diagonalizing the hcrmitian matrix: 



(4) 



where 6b G [— 7r,7r] is the "Blocli angle", and the eigen- 
values Ei{9B) form the band structure of the unper- 
turbed system. The states 



ind£ 



(5) 



where Z^^^ denotes the j-th column of the matrix Z{6b), 
form a complete set of eigenstates of Hq. 

The Landauer conductance g{E) of the disordered wire 
can be calculated by means of the M0ller operator 



n+{E) = 1 



1 



(H - E 



-W G^{E)W, (6) 
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where G^{E) is the retarded Green's function of the dis- 
ordered system. 57+ (E') is an cxd x c» matrix, but its 
non-trivial part has finite dimension. 

If we indicate with Pb the projector on the region of 
length (£ -|- 2) containing the disordered region plus the 
two unperturbed columns located at n = and at n = 
C + 1 and by Pb the projector on the column of index 
£ + 1, the transmission amplitude at the energy E from 
the energy channel 77 on the left, to the energy channel v 
on the right on the disordered region is given by 



A^E,ij^i^) = {E,i^\Pbn+PB\E,y), 



(7) 



where \E, ri{v)) are the Bloch states of energy E belong- 
ing to the ri{v)-ih band (see Eq. ©) and with positive 
group velocity. The conductance is calculated by means 
of the following equation: 

g{E) - 2e^i?(ii;,77)|A-(r; ^ v)\''v{E,u), (8) 



where v{E, v) 



i«i^andi.(iJ,,) = ^(»)' 

are respectively the group velocity on the band v and 
the density of states per column of the underlying per- 
fect wire at the energy E on the band 7?; e is the electron 
charge. We choose the units in such a way that the con- 
ductance g{E) of the perfect wire assumes integer values. 

If the matrix Z is independent of the Bloch angle 9b 
[l^ . and at the energy E all the bands of the perfect 
wire are present, it can be shown that the expression ^ 
assumes the following particularly simple form Q : 







( 
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The meaning of the symbols entering in Eq. jnj is spec- 
ified below: 

• is the standard vector of dimension N having 
one at the i/-th position and zero elsewhere. 



where are the N x N diagonal matrices defined 
by {F^)n^n = e"*^'"^" . Here 0^ are the Bloch angles 
at which the ry-th band has energy E and positive 
(+) or negative (-) group velocity. 

Eq. is the main scattering theory tool for the evalua- 
tion of disorder effects on the conductance of a periodic 
system composed by an arbitrary number N of chains. 
To exemplify its use, we first treat the case of a strictly 
one-dimensional system; the case of two coupled chains 
will be discussed successively. 

For the strictly one-dimensional model {i]=v=\) 
with pure diagonal disorder in the region 
the expansion of the fraction in Eq. ((UJ gives: 



f ^^77 



iXn 



where / = ~i2TTD{E, l)wii and the sums are over all 
the sites within the disordered region; moreover, IS 
equal to if m > n and 9^ ii m < n. Using Eq. © 
the average conductance is 



{g{E))^{\A^\') 



(10) 



Every term contributing to the average conductance can 
be represented by a Feynman diagram. We have found 
that it is convenient to classify all these diagrams into 
two categories: those for which their phase depends from 
the particular sites where the scatterings took place and 
the others. At the lowest order (fourth), there is just one 
non-vanishing diagram belonging to the latter category 
(see Fig. This diagram, in fact, is proportional to the 




m 



• Z, G are block-diagonal matrices with blocks Z and 
G = Z'^GZ, respectively. G is the free Green's 
function projected on the left and on the right of the 
same column. Under the hypothesis just assumed, 
G is diagonal with values Gr,,n = —i2TrD{E, rj). 
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FIG. 1: First diagram contributing to the mean squared con- 
ductance in tlie regime of resonance {6^ — 9~ = tt). 

squares of the disorder parameters (xn,n)^ and (xm,m)^ 
so that, taking the average value it does not vanish; more- 
over, moving either the position n or the position m by 
one or more lattice parameters in either directions, it 
acquires a total phase change equal to e'^^^ ~^ \ The 
mechanism by which the resonance appears as a sharp 
peak is that only for (0"*" — 9~) = qir, with q integer, 
all the amplitudes of such diagrams sum coherently. We 
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name the occurrence of this condition as "tt- couphng of 
the energy bands" . 

It is possible to see that, as far as the condition 
(0+ — 0~) = TT is satisfied, not only the fourth order 
of the average conductance acquires an additional term, 
but also every even-order term of the expansion in the 
disorder strength, so determining the magnification of 
the average conductance. We define the sum of all these 
diagrams as the "resonant contribution to the average 
conductance" {g)res- 

For one-dimcnsional systems the resonance occurs at 
E = 0. With Ti.i = -i (which implies, at E ^ 0, 
Gi_i = —i and 9^ = —6^ = j) and for diagonal disorder 
of intensity wW ~ Wd the sum of all the diagrams up to 
the twelfth order gives: 



{giE = 0)), 
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With similar manipulations for real off-diagonal disorder 
of intensity i ^ Wf wc obtain: 



{9{E = 0)) 



32 

10698752 



1664 
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6! 
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Interestingly, we have noticed numerically that for any 
given value of Wd there exists a value of Wf such that 
the resulting conductances g{E) have almost the same 
off-resonance shape, while the height of the peak for off- 
diagonal disorder is much bigger than the other. In ad- 
dition, we note that for off-diagonal disorder all the di- 
agrams contributing to a fixed order in the expansion 
have the same sign, while for diagonal disorder there is a 
mutual cancellation of diagrams. 

We have extended the previous approach to quan- 
tum wires composed by N different chains all with equal 
nearest-neighbor interactions and containing a disordered 
region always of the same kind and length. We have ob- 
served a strong dependence of the average conductance 
at = from the number N of chains. In particular, 
for N odd and a wide range of real off-diagonal disorder 
0, {g{E = 0; N)) is considerably greater than for the 
even case + 1. By comparing this result with the con- 
ductance evaluated at slightly different energies (where 
the resonance is destroyed), we can conclude that this 
even-odd effect of the conductance is just caused by the 
resonant terms discussed above. 

In the case of more than one coupled quantum chains, 
the rule for calculating the diagrams which are relevant is 
very similar to the one depicted for the one dimensional 
case. The only difference consists in the fact that now 
there are more angles which can propagate and additional 
care has to be taken in selecting the paths contributing 
to the resonance. We give in Eq. H13|) below the fourth 



order expansion of the resonant contribution to the con- 
ductance of the two coupled chains with two resonating 
bands: 



{9)re 



4f(Tr (PiA/„PiM„P2M„P2MnPi 

+P^M^P2Mr,P2M„,P^M^P^ (13) 

+P2M„P2M„PiM,„PiM„P2 
+P2M„PiM„P2M^PiM„P2 

+P2MmPlMnPlMmP2MnP2)) + C.C. 

where n and m arc two arbitrarily chosen different col- 
umn indexes and: 

M, ^ Z^GWu,tZ + P-Z^GWI^Z + Z^GWT,iZF+; 



Pi 



1 




P2 
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The result of the previous trace is in general quite cum- 
bersome, so we report it only for the case where the Bloch 
angles are 9^ = —9^ = ^, and 0^ = —92 — and the 
matrices Z and G are given by: 



cos a — sm a 
sin a cos a 



G = 



-7S 

-73 



(14) 



J we recover the conduc- 



It is easy to verify that for a 
tance at = of two coupled chains with first neighbor 
hoppings all equal to —1. 

Introducing the disorder, of equal intensity, only on the 
hoppings corresponding to the matrix elements (W^)i,i 
and (Wj)2,2, we obtain: 

This sin** a behavior is in good qualitative agreement 
with the numerical data reported in Fig. |21 The most 
interesting cases are those for which the 7r-coupling be- 
tween the bands may occur for energies different from the 
center of the energy spectrum or only between a subgroup 
of all the bands. In this cases, in fact, the (chiral) sub- 
lattice symmetry claimed in the references [7| is broken. 
An example of this occurs for 



U = 



-1 
-1 



T = 



-1 
-0.3 



-0.3 

-1 



(16) 



In Fig. 121 (a) we report the density of states and in Fig.|2| 

(b) the average conductance for the two-chains system 
with off-diagonal disorder. The energy band structure of 
the perfect underlying system is reported in Fig. 01 (c). 
As a further example let us focus on the peak at i? = 1 
where all the bands arc present but only one resonates. 
In this case, an equation similar to Eq. 1131 gives: 

{g)res (X ((«i)2) + {iWl2)') ~ ((^^1^2)'))' + ' ' ' 



4 




FIG. 2: Average conductance for two coupled chains as func- 
tion of the Bloch angle 62 ~ —O2 and of the a angle defining 
the matrix Z in (I14II . The angle 6^ = —9^ are fixed to be 
The peaks appearing for 82=^ are due to the resonance 
contribution and, as function of a, it follows qualitatively the 
sin* a predicted by 11511 . A gray-scale has been used to em- 
phasize the peaks. 



numerical simulations [g]. 

From the evaluation of the spectrum of eigenvalues for 
several configurations we have obtained the average DOS 
depicted in Fig. |21(a). This figure evidences a direct con- 
nection between the peaks in the conductance and in the 
density of states. As far as wc know, such a connection 
has been demonstrated for strictly one dimensional sys- 
tems only 0. 

Finally we justify why we kept in the development of 
the mean conductance only terms of the form w^'^ with 
p = q and no smaller values of q. The reason is that 
the average conductance {g{E] w, C)) converges at a fixed 
non-trivial value for every energy when C 00 with w^C 
constant {w represents the global strength of disorder). 
In practice we have verified that this convergence is very 
fast Hj. 

We would like to thank A. Cresti for having provided a 
very fast code, based on the Keldysh formalism i used 
as a check of our code and as a source of large statistics 
of the conductance. 



LU 
_1 

O 
z 

< 
X 

o 
o 

_1 

CD _ 



I \ 


(b) \^ 




(c) 





0.6 



0.4 r 



-5^-3-2-1 1 2 3 4 
ENERGY 



FIG. 3: Average DOS in arbitrary units (a), average conduc- 
tance (b), and band structure(c) of the two coupled chains 
defined by . The density of states corresponds to the dis- 
ordered region (device) isolated from the leads. There are 
three energies at which 7r-coupling occurs. In correspondence 
of these energies, the average conductance and the density of 
states are peaked. 



The latter expression shows clearly that these two kinds 
of off-diagonal disorder (on the T and on the U matrices) 
compete with each other, and this is consistent with our 
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